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ABSTRACT 

Massive early- type galaxies are observed to lie on the Mass Plane (MP), a two- 
dimensional manifold in the space of effective radius i? c , projected mass M^ 2 (measured 
via strong gravitational lensing) and projected stellar velocity dispersion <t o2 within 
R e /2. The MP is less 'tilted' than the traditional Fundamental Plane, and the two have 
comparable associated scatter. This means that the dimensionless structure parameter 
c o2 = 2GM^ 2 / (R c (j^ 2 ) is a nearly universal constant in the range cr c2 = 175 — 400 km 
s . This finding can be used to constrain the mass distribution and internal dynamics 
of early-type galaxies: in particular, we explore the dependence of c o2 on light profile, 
dark-matter distribution, and orbital anisotropy for several families of spherical galaxy 
models. We find that a relatively wide class of models has values of c e2 in the observed 
range, because c e2 is not very strongly sensitive to the mass distribution and orbital 
anisotropy. The degree of fine tuning required to match the small intrinsic scatter of c c2 
depends on the considered family of models: if the total mass distribution is isother- 
mal (cx r~ 2 ), a broad range of stellar luminosity profile and anisotropy is consistent 
with the observations, while Navarro, Frenk & White dark-matter halos require more 
fine tuning of the stellar mass fraction, luminosity profile and anisotropy. If future 
data can cover a broader range of masses, the MP could be seen to be tilted by the 
known non-homology of the luminosity profiles of early-type galaxies, and the value of 
any such tilt would provide a discriminant between models for the total mass-density 
profile of the galaxies. 

Key words: galaxies: elliptical and lenticular, cD - galaxies: formation - galaxies: 
kinematics and dynamics - galaxies: structure - gravitational lensing 



1 INTRODUCTION 

The origin of empirical scaling laws is a key open issue 
in observational cosmology. Galaxies do not come in all 
sizes, shapes, colours, but rather tend to to live in lower- 
dimensional manifolds, which represent a stringent testing 
ground for theories of galaxy formation and evolution. 

Early-type galaxies obey a particularly tight 
scaling law: the so-call e d Fundament a l Pl ane (FP; 
iDiorgovskv fc Davis! 1 19871 : iDressler etabl Il987f ). In the 
space of effective radius R c , central velocity dispersion 
<T e 8 and effective surface brightness (I) c =L/(2tvR 2 ) (where 
L is the total luminosity of the galaxy), they lie on the 
foll owing relation with re markably small scatter ( < 20% in 
R.: lBernardi et al.|[2003bl ): 



where the numerical value of a and b depends somewhat 
upon the wavelength of o bservations and upon the sample 
and t he fitting method |Pahre. Diorgovski fc de Carvalhol 
Il998l ; iBernardi et al.l l 2003bl) . The FP is said to be 'tilted', 
in the sense that the coefficients a and 6 differ signif- 
icantly from the values a — 2 and b = —1 expected 
for structurally and dynamically homologous systems with 
luminosity-independent stellar mass-to-light ratio and dark- 
matter distribution. Several explanations have been pro- 
posed for the tilt, including a systematic dependence of 
stellar mass-to-light ratio or dark-matter content and dis- 
tribution upon luminosity (and hence presumably upon 
mass) , structural non- h omology and orbital anisotropy 



log 7i e = alog<j c8 + 61og(J) c + const, 



(1) 



(e.g. iFaber et al.1 Il987l ; iBender. Burstein fc Faberl 
Renzini fc Ciottil 1 19931: ICiotti. Lanzoni fc Renzinil 



Borriello. Salucci fc Danesel feooa) 
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Recently, iBolton et all (|2007l . l2008bl ). using a sample 
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of strong gravitational lenses, have shown that early-type 
galaxies lie on a Mass Plane (MP) 



log R c = a m log (j C 2 + bm log E e 2 + const , 



(2) 



where er 2 is the projected velocity dispersion within an aper- 
ture radius R e /2 and E C 2 is the surface mass density within 
R c /2, with a m = 1.82 ± 0.19, b m = -1.20 ± 0.12 and RMS 
orthogonal scatter of 1.24 when normalized by the observa- 
tional errors. The fact that (a m , 6 m ) are close to (2, —1) and 
that the scatter is small can be expressed in terms of struc- 
tural and dynamical homology of the lenses, by defining the 
dimensionless structure parameter 



Ce2 



2GM£ 



(3) 



where M% 2 is the total projected mass within Re/2. For their 
sample of lens ea rly- type galaxies from the Sloan Lens ACS 
(SLACS) Survey iBolton et al.1 (|2008bl ) find on average 



(logc e2 ) = 0.53 ±0.057, 



(4) 



which throughout the paper we will refer to as the "observed 
range" of c C 2. We note that the observed scatter on (log c e 2) 
is 0.08, b ut here we consider t he estimated intrinsic scatter 
0.057 (see IBolton et al.ll2008bl ). 

As discussed in several papers (IBolton et al. 



2006; 



iTreu et al.ll2006l ; IBolton et ai1l2008al ; iTreu et al.ll2008l ) the 

SLACS lenses are found to be indistinguishable from con- 
trol samples of Sloan Digital Sky Survey (SDSS) galaxies 
with the same stellar velocity dispersion and size, in terms 
of luminosity /surface brightness, location on the FP, and 
environment. This inspires some confidence that the results 
found for the lens sample, including the MP, are generic 
properties of the overall class of early-type galaxies. 

Independent of its origin and theoretical interpretation, 
the existence of the MP is a powerful empirical tool to esti- 
mate galaxy ma ss by using informa tion on size and velocity 
dispersion only (|Bolton et al.l [20071 ). In addition, it is clear 
that the very existence of the MP may be used to improve 
our understanding of galaxy formation and evolution. What 
is the origin of such a strong correlation among measurable 
galaxy quantities? Or, in other words, what kinds of galaxy 
models can be ruled out by the existence of a tight MP? 
Although this question has been asked before in regards to 
the traditional FP, the MP provides an additional power- 
ful tool. In fact, there are a few differences between the FP 
(equation [1} and the MP (equation [2]) : 

(i) The FP is sensitive to the galaxy stellar mass-to-light 
ratio, while the MP is not. This implies that, e.g., the role 
of stellar populations in establishing the tilt and scatter of 
the FP can be disentangled by looking at the MP0 

(ii) The FP is traditionally based on the central velocity 
dispersion a c s (measured within R c /8), while the MP has 



1 Strictly speaking the MP depends on the properties of the stel- 
lar populations through _R C and a e 2 . However, R e and cr c 2 do not 
depend on the value of the stellar mass-to-light ratio, but only on 
its radial variation. This variation is expected to be small based on 
observed col our gradients and is generally neglected in dyna mical 
studies fe.g. lKronawitter et al"1l2000l ; ICappellari et al.ll2006l ). For 
simplicity, in this study we assume uniform stellar mass-to-light 
ratios within each galaxy. 



been constructed using <r C 2, which is measured within R e /2. 
This is a consequence of the fixed spatial observing aperture 
of the SDSS spectrograph; an MP based upon a a g could 
be constructed using spatially resolved spectroscopy of the 
SLACS lens sampleQ 

(iii) The FP combines quantities evaluated on different 
scales (_R e , R e /8), while MP combines quantities evaluated 
within the same radius R c /2. Again, this is partially due to 
the fixed SDSS spectroscopic aperture, though the apertures 
of the lensing mass measurements are fixed by the cosmic 
configuration of the individual strong-lens systems. 

Each of the points above can contribute to make the MP 
less tilted (and presumably with less scatter) than the FP. 
For example, given the relatively large spectroscopic aper- 
ture used to define c C 2, we expect it to be robust with re- 
spect to changes in the detailed properties of galaxy struc- 
ture, internal dynamics, and dark- matter content. Similarly, 
by replacing surface brightness with surface mass density 
we expect that tilt and scatter due to diversity of chemi- 
cal composition or star formation history be reduced in MP. 
Furthermore, having all but removed the effects of stellar 
population the MP is potentially a cleaner diagnostic than 
the FP of the structural and dynamical properties of early- 
type galaxies. 

In this paper, we exploit the existence of the MP 
to constrain important properties of early-type galaxies, 
such as orbital anisotropy and dark-matter distribution. 
We achieve this goal by constructing observationally and 
cosmologically motivated families of galaxy models and 
finding the range of parameter spaces consistent with 
the observed range of c C 2. For the sake of simplicity, in 



symmetric models. As with the FP dFaber et al 



Sagl ia. Bender fc Dressle"rl 19931; Prugniel fc Simien 





1987 




1994 



spherical symmetry is expected to increase the scatter of the 
MP because of projection effects. Thus, a natural follow-up 
of the present work would be the extension to non-spherical 
models. 

The paper is organized as follows. In Section [2] we de- 
scribe the models, in Sections[3]and|4]we present our results, 
and in Section [5] we conclude. 



2 MODELS 

2.1 Methodology and general properties 

We consider spherical galaxy models with stellar density 
distribution p*(r) and total density distribution ptot(f). 
The radial component ay(r) of the velocity dispersion 
tensor is given by so lving the Jeans equation (e.g. 
iBinnev fc TremaindEfflOsI ) 



dp t a1 2 / 3p»o- 2 



dr 



+ 



= ~P*9r, 



(5) 



2 In general, the larger the aperture radius R used to measure 
the aperture velocity dispersion tr 2 , the less ct 2 is sensitive to 
the orbital anisotropy. We recall that for any stationary, non- 
rotating, spherically symmetric system with constant mass-to- 
light ratio cr 2 (R) — > o \ 2 /3 for R — > oo, where <r 2 is the virial 
velocity dispersion (e.g. ICiottilll994h . 
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where g r (r) = d$(r)/dr, $(r) is the total gravitational po- 
tential generated by ptot(r), and 



14 



P{r) 



2a 2 . 



(6) 



is the anisotropy parameter (a$ and o"^ are, respectively, the 
# and <y3 components of the velocity-dispersion tensor). 

The line-of-sigh t velocity dispersion is 

(jBinnev fc Mamon|[l983 ) 



a 2 p (R) = 



where 



MR) 



d2 



p*(r)a 2 rdr 
Vr 2 - R 2 ' 



p*(r)rdr 
Vr 2 - R 2 



(7) 



(8) 



is the stellar surface density (we assume that the stellar 
mass-to- light ratio is independent of radius). The aperture 
velocity dispersion within a projected radius R, the closest 
analog to the measured stellar velocity dispersion, is deter- 
mined via 



l(R) 



where 



2tt 



M*(R) 



/ T, t {R')al(R')R' dR', 
Jo 



Mf(R) = 2tt MR')R' dR' 



is the projected stellar mass within R. So 



(9) 



(10) 



l(R e /8) 



and <T c2 = <T a (i? c /2). Note that the mass weighting expressed 
here is equivalent to luminosity weighting for the case of a 
spatially uniform stellar mass-to-light ratio. 

Gravitational lensing analysis allows one to measure the 
total projected mass density within the Einstein radius. The 
total projected mass within a radius R of a spherical system 



M t p ot OR) = 27T f Etot (R')R' dR', 
Jo 



where 



Etot (B) = 2 



Ptot{r)rdr 



R 2 



(11) 



(12) 



is the total surface density. The mass within the Einstein ra- 
dius, measured by gravitational lensing is obtained by set- 
ting R equal to the Einstein radius. The size of the Ein- 
stein radius depends on the geometry of the lensing sys- 
tem, through the angular diameter distances between the 
observer, lens and background source, as well as on the mass 
distribution of the lens. Typically, for galaxy-size lenses, Ein- 
stein radii are of order of one arcsecond, or ~5 kpc for lenses 
at moderate redshift. For the SLACS lens sample, the Ein- 
stein radii are typically about half the effective radius of the 
lens galaxy. 

For a given projected radius R we define the structure 
parameter 



c(R) 



GMf ot {R) 
Rai(R) ' 



(13) 



So, by definition (equation [3} c C 2 
ACt(i?c/2). 



c(R e /2), because M| 
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Figure 1. The ratio £ = r B /R c as a function of the effective radius 
R e according to equation 1124 \ for NFW plus stars and adiabati- 
cally contracted NFW plus stars models (see text for details). 

2.2 Stellar density distribution 

We consider two families of stellar density distributions — 7 
models and Sersic models — that are known to match well 
the observed surface brightness profiles of early-type galax- 
ies over the range of interest ~1-10 kpc for constant stel- 
lar mass-to- light ratios. The density profile of the 7-models 
jDehnenll 19951 ; iTremaine et al.lll994 ) is given by 

P*ir)= \J ^ r „ M * r * (0< 7 <3), 



4tt r~t(r t +r) 4 



(14) 



where M» is t he total ste l lar m ass. The cases 7 = 1 and 
7 = 2 are the iHernquistJ (|l990h and Ijaffd (|l983h models, 
respectively. 

In the case of the Sersic models , the projected density 
distribution follows the lSersicI (|l968t l R 1/m law: 



E.(.R) = E.,o exp 



^ (I 



l/m' 



(15) 



where b(m) ~ 2m- 1/3 + 4/ (405m) (jCiotti fc Bertinlll999l ). 
For m = 4 the de Vaucouleurs' (1948) R 1 ^ 4 law is obtained. 
By deprojecting E* one obtains the correspond ing intrinsic 
density distribution (jBinnev fc TremaindlBoSn 



7T 



dE, dR 



dR y/R? 



(16) 



2.3 Total density distribution 

We consider four different models for the total den- 
sity distribution. First we consider a singular isother- 
mal sphere (SIS) model, which provides a generally 
good descripti on of the lensin g properties of early-type 



galaxies (e.g. Kochancj Il994j : iTreu fc Koopmansl 120041 ; 



iKoopmans et al.1 l2006h. as well as of their stellar kine- 
matics (e.g. iBertin etal] Il994l ; iGerhard et~ail l200lh . Sec- 
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Figure 2. Total (solid lines), stellar (dotted lines), and dark- 
matter (dashed lines) projected mass distributions M P (R) for 
m = 4 Sersic models with non-contracted (top) and adiabatically 
contracted (bottom) NFW dark- matter halos of given § = r E /R e 
and stellar-to-total mass fraction /». The projected mass M P (R) 
is in units of M, +M(j m . All models have concentration C200 = 8. 



ond, we consider light-traces-mass (LTM) models. Although 
they are known to fail observational constr aints based 
on st atistical analyses of stron g lenses (e.g. , Rusin et alj 
120031 : iRusin fc Kochanekl I2005I: iBolton et all l2008bl ). stel- 
lar dynamics (e.g., iGerhard et alj|200ll). individual strong- 
lens mode l ing (e . g., IWavth et all 2004 iDve fc WarreifeOOa 
I Dye et all 120071: Gavazzi et al.l 20081 ) , weak-lensing anal- 
ysis (e.j 



, Gavazzi et al.l 20071 ) , and combined strong- 
lens/dynamical analyses (e.g., iKoopmans fc Treul |2002| . 



man" hypothe- 
we consider two 



2003) : iTreu fc Koopmandliooil , 120031 , 12004 IKoopmans et al.1 
2006), they serve as a useful "straw 
sis to test against the MP. Lastly, 
families of cosmologica ll y mo tivated models based on 
iNavarro. Frenk fc White] |l996l , NFW) halos with the ad- 
dition of stars. In one case (NFW plus stars models, here- 
after NFW+S) the stars are added leaving the halo unper- 
turbed; in the other one (adiabatically contracted NFW plus 
stars models, hereafter acNFW+S), the halo is assumed 
to "respond" to the sinking of baryons towards the cen- 
tre of the galaxy as prescribed by the adiabatic contraction 
recipe of iBlumenthal et al.l (Il986l). Sev eral arguments sug- 



gest that the iBlumenthal et aL (1980} mode l might over- 
estimate the compression of the halo (se e iGnedin et al.l 
2004 lEl-Zant et al.ll200i : iNipoti et al.ll2004l ). In particular, 



Gnedin et al.l (|2004) argued that the s tanda rd adiabatic con- 
traction model of IBlumenthal et al.l (|l986l ), based on some 
simplifying assumptions such as spherical symmetry and cir- 
cularity of the particle orbits, tends to overpredict the in- 
crease of dark matter density in the central regions. How- 
ever, considering both NFW+S and acNFW+S we should 
bracket the realistic ran ge of NFW halo models (see also 
I Jiang fc Kochanekll2007l ). 

We now define for each model the total (stellar plus 
dark matter) density profile ptot(r). As we limit to spher- 
ically symmetric distributions, in all cases the modulus of 
the gravitational field is 



9r(r) = 



GM tot (r) 



(17) 



where M t ot(r) = 4-7T J Q ptot(r')r' 2 dr' the total mass within 
r. 

In the case of SIS models, the total (stars plus dark 
matter) density profile is 



ptot(r) = 



(18) 



2nGr 2 ' 

where asis is the one-dimensional velocity dispersion of 
isotropic SIS. For LTM models the total density profile is 
Ptot(?") = Kp*( r ), where K is a dimensionless constant (note 
that c 2 is independent of K). In NFW+S models the dark- 
matter halo is described by a NFW model, so the dark- 
matter density distribution is 

M dm 1 



Pdm(f) = 



(19) 



47r/(C 2 oo) r(r + r s ) 2 ' 

where r s is the scale radius, and the distribution is trun- 
cated at the virial radius rvi r . The average dark-matter den- 
sity within 7Vi r equals 200 times the critical density of the 
Universe. In the equation above, C200 = Pvir/r B is the con- 
centration parameter, 

C200 



/(C200) =ln(l + C 20 o)- 



(20) 



1 + C200 

and Afdm is the total dark-matter mass. In this case the total 
(stars plus dark matter) density profile is 

Ptot(r) = pdm(r) + p,(r). (21) 

Formally, these NFW+S models have three free parame- 
ters: concentration C200, the stellar mass fraction /» = 
M*/(M« + Mdm) and the ratio £ = r s /R a . Cosmological 
simulations sugge st values of C200 ~ 7 — 9 for low-redshift 
galaxy-size halos (|Neto et al.ll2007l ). Thus, we fix C200 = 8, 
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and we explore different combinations of values of £ and /„ . 
There are indications th at f* is typically of few p ercent in 
early- type galaxies (e.g. I Jiang fc Kochanekl 120071 ). so here 
we consider the cases /* = 0.02 and /» = 0.1. The value of 
£ for given /» is not strongly constrained by models and ob- 
servations, but for the present investigation it is sufficient to 
individuate a realistic range of values of £. For this purpose 
we can use the observed correlation between effective radius 
and total stellar mass for early-type galaxies: 



Rc 
kpc 



= B 



So 



where A = 0.56 and B = 2.88 x 10" 
By definition of the virial radius, 

3M dn 



47T7-3 



200p c 



(22) 

(|Shen et al. Il2003h . 

(23) 



where p cr i t — SHq / (8nG) is the critical density of the 
Universe and Ho is the Hubble constant (here we assume 
Ho = 70kms~ Mpc^ 1 ). Combining equations (|22[) an d 
()23p . we get the following relation between £ and R e 



1 - /. R 



(l/A)-3 



8007rSV^ f„Cl 



Pcrit 



(24) 



where R c = i? e /kpc and p crit = p CI it/(M Q kpc~ s ): £ as 
a function of R e is plotted in Fig. [T] Early-type galax- 
ies of the SLACS sam ple have effective ra dii in the range 
0.3 < log R c / kpc < 1.3 (|Bolton et al.ll2008aT ). For each value 
of /„ we show results for two values of £ roughly correspond- 
ing to the upper and lower limits of this range. In particular, 
we consider £ = 6 and £ = 13 when /» = 0.02, and £ = 3 
and £ = 8 when /» = 0.1. The resulting stellar, dark-matter 
and total projected mass profiles are plotted in Fig. [5] (top) 
when the stellar profile is a de Vaucouleurs (or m — 4 Sersic) 
model. 

We also consider acNFW+S models, in which the dark- 
matter halo is adiabatically contracted following t he stan- 
dard recipe (|Blumenthal et al.|[l98r3 ; iKeetonl l200lh . In the 
considered models, the dark-matter profile is obtained by 
adiabatically compressing an initial NFW profile with the 
same values of the parameters C200, £ and /„ as for the 
non- compressed models. The final dark-matter distribution 
is computed numerically for each given stellar distribution. 
Projected mass profiles of acNFW+S models with de Vau- 
couleurs' stellar density distribution are shown in Fig. [5] 
(bottom) to allow a direct comparison with the correspond- 
ing non contracted models shown in the top panels. 



2.4 Orbital anisotropy 

We consider two parameterizations of radial anisotropy in 
the stellar di stribution: constant anisotro py and Osipkov- 
Merritt COM; IOsipkovlll979l ; lMerritdll985( ). In the first case 
the value of the anisotropy parameter is the same at all radii: 



const, 



(25) 



an d the radial component of the velocity dispersion tensor 
is (jBinnev fc Tremaindl2008h 



2 ( > 

a r (r) = 



,.2ii 



p*{r) 



„/20„ r^dHr') , 
dr> 



P*( r )- 



(26) 



In the case of OM anisotropic models, the anisotropy in 
the stellar orbital distribution is introduced by using the 
following parameterization: the radial dependence of the 
anisotropy parameter is 



r 2 + ri 



(27) 



where the quantity r a is the so-called "anisotropy ra- 
dius". For r 3> r a the velocity dispersion tensor is radially 
anisotropic, while for r<r a the tensor is nearly isotropic. 
Isotropy is realized at the model centre, independently of the 
value of r a . In the case of OM models, the radial component 
of the velocity dispersion tensor is given by 



p*(r)o-r(r) = 



ri + r 2 



where 



Q*{r) 



1 + 



p*(r) 



(28) 



(29) 



(|Merrittlll985h . 

We consider constant-anisotropy models because they 
are the simplest possible anisotropic models, and span 
the full range of anisotropies, from tangentially to radi- 
ally biased. However, OM models should be more real- 
istic, because observational indications suggest that typi- 
cal massive elliptical galaxies are, in the c entral regions, 
isotropic or mildly radially anisotropic (e.g. I Gerhard et alj 
l200ll : ICappellari et al.1l2007h . and different theoretical mod- 
els of galaxy formation predict that elliptical galaxies 
should have anisotropy varying with radius, from almost 
isotropic i n the centre to radially bia s ed in the out- 
skirts (e.g. Ivan Albadalll982l: lBarnesill992l : lHernquist|[l993l : 
iNipoti. Londrillo fc Ciottill2006h . 



2.5 Consistency and stability 

A galaxy model is consistent if it has a positive distribu- 
tion function. Not all combinations of the parameters in- 
troduced in the sections above generate consistent models. 
For instance, a necessary (bu t not sufficient) conditio n for 
consistency of OM models is (jCiotti fc Pellegrimlll992l ) 



< 0, 



(30) 



dg*(r) 
dr 

which holds also in the presence of a dark-matter halo. For 
{3 = const models, the necessary con dition is 7* > 2/3 , 
where 7* = — lim r _>o (dlnp*(r)/dlnr) l| An fc Evansl 120061 . 
see also Richstone & Tremaine 1984 and Tremaine et al. 
1994). This condition holds not only for one-component sys- 
tems, but also for two-co mponent systems if p to t ~ r~' 2 ~ E ' 
for r -> 0, with e > (|An fc Evansl 120061 ). In the mod- 
els here considered e > 0, except for the limiting cases of 
7 = 2 models and/or SIS total density, in which e — 0. How- 
ever, lAn fc Evansl (|2006l) show that the necessary condition 
is 7* > 2/3 + (i - /3)n if ptot oc r- {2+n) , with n > 0, so 
one can argue that the condition 7* > 2/3 must hold also if 
ptot oc r~ 2 for r — > (i.e., logarithmically divergent central 
potential). 

Thus, the requirement of consistency reduces the pa- 
rameter space. To ensure physically meaningful models we 
computationally check consistency and rule out regions of 
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parameter space that would give rise to non-consistent mod- 
els. In particular, we exclude OM models that do not satisfy 
the condition (|30[) . and f3 — const models with 7* < 2/3. We 
recal l here that for Sersic models 7, = (m — l)/m (|Ciottil 
1991), while for 7 models simply 7* = 7. 

Additional constraints would come from the require- 
ment of model stability. In particular, strongly radially 
anis otropic systems are expected t o be radial-orbit unsta- 
ble (|Fridman fc Polvachenkol Il984f ). However, while there 
are robust estimates of the maximum amount of radial or- 
bital anis otropy allowed for stable one-component systems 
(see, e.g.. iMerritt fc Aguilarlll985l; iBertin fc Stiavellij Il989l ; 



ISahal Il99ll : iMeza fc Zamoranol 1 19971 ). 

much less is known 

about the stability of two-component systems, though there 
are indications that the presence of a massive halo con- 
tribu t es to the systems' stability (e. g. IStiavelli fc Sparkel 
Il99ll ; iNipot i. Londrillo fc Ciottil |2003 ). As a consequence, 
we have not enough information to exclude models on the 
basis of stability arguments^ Nevertheless, one needs to 
bear in mind that models with extreme radial anisotropy, 
even satisfying the necessary consistency condition, might 
be non-consistent or radially unstable. 



3 RESULTS: CONSTRAINING MODEL 
PARAMETERS WITH OBSERVATIONS 

We numerically compute c C 2 for 7 models with < 7 < 2 
and Sersic models with 1 < m < 10 for different total 
mass distributions and orbital anisotropy. Here we discuss 
how the obtained values of c C 2 compare with the observed 
range (logc e2 ) = 0.53±0.057 (jBolton et alj|2008bl ). For each 
combination of family of stellar systems and total density 
distribution, we represent our results using contour plots 
of log c C 2 in planes of anisotropy versus stellar profile pa- 
rameters (Figs. I3I8|I : /3-7 for /3 = const 7 models, f3-m for 
/3 — const Sersic models, log(r a /i? e )-7 for OM 7 models, and 
log(r a /i? e )-m for OM Sersic models. In such plots, "allowed" 
regions are white, regions corresponding to non-consistent 
models are dark-shaded, while regions outside the observed 
range are light-shaded. 

OM models that are unphysical because non-consistent 
(dark-shaded areas in the right-hand columns of Figs. [3} 
[5]) are also outside the observed range. On the other hand, 
there are non-consistent ft — const models with c C 2 in the 
observed range. This different behaviour of OM and (3 = 
const models is not surprising, because the latter are known 
to be less realistic than OM models (see Section[2]4]|, and, in 
particular, some radially anisotropic /3 — const models turn 
out to be unphysical because they are radially anisotropic 
down to the very centre of the system. This finding stresses 
the importance of investigating consistency when modeling 
observational data. 

At a first level of interpretation, the plots in Figs. 1 3181 
show that a relatively wide class of models have values of the 
structure parameter c C 2 in the observed range. Thus, the fact 
that the observed values of c e 2 lie in a small range does not 



3 All the models considered in this paper are two-component 
models. In LTM models we assume that the dark matter has 
the same density distribution as the baryons, but not necessarily 
the same velocity distribution. 



necessarily imply that early-type galaxies are structurally 
and dynamically homologous. The reason for this is that c e 2 
is not dramatically sensitive to the mass distribution and 
orbital anisotropy. However, a more detailed analysis of the 
diagrams indicates that there is also a wide class of mod- 
els that lie outside the region allowed by the observations. 
Here we summarize the behaviour of families of models with 
different total mass distribution: 

• SIS models (Fig. [3j|: SIS models are consistent with 
observational constraints for a wide class of stellar density 
profile and anisotropy. All isotropic 7 models and isotropic 
Sersic models with m<5.5 have c C 2 within the observed 
range. Higher-m Sersic models can be reconciled with the 
observations if their stars have radially-biased orbits. Strong 
radial and tangential anisotropy is excluded. However, very 
strong radial anisotropy should be excluded on the basis of 
consistency arguments, while, as briefly discussed in Sec- 
tion [23] very strong tangential anisotropy is not expected. 
One might speculate that the OM models excluded by the 
observed range trace roughly the region of radially unstable 
models (see Section [2.51) . 

• LTM models (Fig. [3J|: isotropic (and mildly radially 
anisotropic) 7 models with LTM potential are consistent 
with the observational constraints, while isotropic (and 
mildly radially anisotropic) Sersic models with LTM poten- 
tial are acceptable only for m < 5.5. In contrast with the case 
of SIS models, there is no way of reconciling m>6.5 Sersic 
models with the observations. As in the case of SIS models, 
the lower limit on r a /i? for OM models might be deter- 
mined by radial-orbit instability. More tangential anisotropy 
than in SIS models is allowed, though with some fine-tuning 
with the stellar profile parameters m and 7. Curiously, if 
one considers 7 = 2 models under the assumption of LTM 
and constant anisotropy, the observational constraints would 
favour tangential with respect to radial anisotropy. 

• NFW+S models (Figs. [5] and 16}: for some NFW+S 
models only remarkably small regions of the parameter 
space are allowed, in contrast with the case of SIS mod- 
els. The worst case is that of more dark-matter dominated 
models (/* = 0.02, £ = 6): remarkably, isotropic 7 mod- 
els are excluded, and isotropic Sersic models are allowed 
only for m > 7.5). Better is the most baryon dominated case 
(/« = 0.1, £ = 8), which — as expected — behaves similarly to 
the LTM case, so OM Sersic models can be reconciled with 
the observations only for m<6. In some cases (e.g. OM 
Sersic NFW+S models) the anisotropy and stellar-profile 
parameters must be fine-tuned in order to have c C 2 in the 
observed range. For NFW+S models with acceptable c e 2 we 
find 0.1 < /dm < 0.6 and 0.2 < f» m < 0.7, where / dm and f* m 
are, respectively, the intrinsic and projected dark-matter-to- 
total mass ratios within R a . 

• acNFW+S models (Figs.[7]and|8]): overall, adiabatically 
contracted models behave similarly to non-contracted mod- 
els, but allowed regions in the parameter space are slightly 
more extended in acNFW+S models than in the correspond- 
ing NFW+S models. As in the case of NFW+S, baryon- 
dominated models (larger values of /» and £) are more suc- 
cessful than dark-matter dominated models (smaller values 
of /* and £). However, acNFW+S have some of the same 
undesirable features present in NFW+S models, such as 
a wide class of unacceptable isotropic models (especially 



The Mass Plane of early-type galaxies 7 



when /, = 0.02). For dark-matter dominated OM Sersic 
acNFW+S models to have c G 2 in the observed range, the 
anisotropy and stellar-profile parameters must be fine-tuned. 
The acNFW+S models with e C 2 within the observed range 
have 0.3 < f dm < 0.65 and 0.35 < /£ m < 0.75. 

Summarizing, our results indicate that the tightness of 
the MP requires some degree of "fine tuning" in the inter- 
nal properties of early-type galaxies. Although no family of 
models is strictly ruled out, some families of models allow 
for more freedom in the remaining parameters describing, 
e.g., the luminous profile and stellar orbits. The tightness 
of the MP is well consistent with the hypothesis that the 
total mass profile is that of a SIS, allowing a broad range of 
values for the other parameters. The LTM hypothesis is ac- 
ceptable if early-type galaxies have intrinsic stellar density 
profiles described by 7 models or Sersic projected stellar den- 
sity profiles with index m < 5.5. For either non-contracted 
or adiabatically contracted NFW cases, baryon-dominated 
models require less fine-tuning than dark-matter-dominated 
models. In general, isotropic or mildly radially anisotropic 
velocity distribution is easier to reconcile with the MP than 
tangential or extremely radial anisotropy. 



4 SERSIC INDEX AND THE TILT OF THE 
MASS PLANE 

In the previous Section we have interpreted the average value 
and intrinsic scatter of the parameter c C 2 describing the MP. 
In this Section we consider how future observations of the 
MP over a broader range of galaxy masses could be used to 
further constrain the internal structure of early-type galax- 
ies, based on the known structural non-homology of their 
luminous component. 

The Sersic index m of early-type galaxies correlates 
with galaxy size, in the sense that more extended galax- 
ies have higher m |Caon. Capaccioli fc D'Onofrid Il993l ; 
ICraham fc Guzmanll2003l; IFerrarese et alj|2006ft ~ in particu- 
lar, ICaon. Capaccioli fc D'Onofrid (|l993ft found 



log m = 0.28 + 0.52 log 



Re \ 

kpc J 



(31) 



As the effective radiu s increases with galaxy luminosity (e.g. 
Bcrn ardi et al ] |2003ah . equation (|31|l should imply a correla- 
tion between m and luminosity. For example, using surface 
photometry from the well-defined Virgo Cluster sample of 
IFerrarese et al.l (|2006ft , we obtain the following relation 



logm = (0.27 ±0.02) log 



Lb 



,0 



(2.07 ±0.02) 



(32) 



where Lb is the B-band lumino sity obtaine d from 
the ACS filters as described by iGallo et al.l <|2008h . 
and the uncertainties on the best-fit coefficients are 
1 a (the RMS scatter about the relation is 0.16 
dex). Such a systematic variation of m with luminos- 
ity could at least partly contribute to the observed 
tilt of the FP of early - type g alaxies dHiorth fc Madsenl 



arly - 

19951; ICiotti fc Lanzonl Il997l; ICraham fc Collessl Il997l; 



Bertin. Ciotti fc Del Principel2002l;lTruiillo, Bur kert fc BeU 



20041 ; Cappellari et al.ll2006n . If the stellar mass contributes 



significantly to the projected mass within R c /2 — as is gen- 
erally assumed and as supported by the modeling results 
above — the dependence of m on luminosity could introduce 
a tilt in the MP. This would in turn induce a tilt in the rela- 
tion between M p 2 = M p t (i?. c /2) and M c d 2 im = i?. c a c 2 2 /(2G), 
because the structure parameter c e 2 depends on m for given 
total mass distribution and orbital an isotropy. For the mas - 
sive galaxies studied in their sample, iBolton et al.l (|2008bh 
found 



log 



M? 2 



lO^M, 



= 8 



10 11 M P 



± log c c2 , 



(33) 



with logc e2 = 0.54±0.02 and 5 = I.03±0.04. Note that here 
we are not considering the average value of log c G 2 as in equa- 
tion Q, but rather are accounting for a possible dependence 
of c C 2 on mass. A deviation of S from unity is a signature 
of tilt, so the ob servational da t a are c onsistent with absence 
of tilt. Actually IBolton et al.l (|2008bh found no correlation 
between galaxy mass/luminosity and Sersic index in their 
sample, consistent with the fact that the SLACS sample is 
confined to a relatively small range of galaxy luminosities 
towards the bright end of the luminosity function of early- 
type galaxies: the mass range of Bolton et al.'s sample is 
10.3 < log(M e p 2 /M ) < 12. We thus compute the tilt to de- 
termine whether it is measurable and could lead to further 
discriminatory power if one considered a sample covering a 
larger mass range. For this purpose, we use our models to 
quantify the tilt introduced in the MP by the expected de- 
pendence of m on mass. We consider our models with Sersic 
stellar distributions and we map m into M by combin- 
ing eq uation (|32p with the best-fit correlation ((Bolton et alj 
2008]]) 



log 



Ly 



10 n L 



V0 



= 0.73 log 



K2 



0.24, 



(34) 



assuming B -V = 0.96 (|Fukuaita et alj 1 1995ft . In order to 
isolate the effect of the structural non-homology of the stel- 
lar distribution, it is useful to compare models under the 
same assumptions on the velocity distribution. In this Sec- 
tion we focus on models with radial anisotropy with OM 
parameterization. We consider two families of OM models 
with values of r a /_R c independent of m: isotropic models 
(r a /i? c = 00) and radially anisotropic models with r a = R c . 
As for the FP, dynamical non-homology might contribute to 
produce a tilt of the MP, so it is important to consider also 
the case in which more massive (higher-m) system are more 
radially anisotropic than less massive (lower-m) systems. 
This choice is motivated by the fact that higher-m system 
can sustain more radial anisotropy than lower-m systems 
(see Sect ion | 2.5land Figs.[3]|8l) and by prev i ous studies on the 

Lanzonil 



Ciotti. Lanzoni fc Renzin 



>y prcvior 
i 1096; C 



iotti 



FP ( 

1 1997i ; iNipoti. Londrillo fc Ciottil I2OO2T ) . Thus, we also con- 
sider a family of models in which log r a /_R e = 0.1(10 — 
m) — 0.3: with this parameterization m ~ 1 systems are 
almost isotropic, while m ~ 10 systems are strongly radially 
anisotropic. In Fig.[9]we plot M p 2 /Af^ 2 lm as a function of A/ p 2 
(and of m) for OM R}'™ models with different total mass 
distribution and for the three different choices of the value 
of r&/Re (isotropic in the left-hand panels, r a = R c in the 
central panels and r a /i? e depending on m in the right-hand 
panels). In each panel the dashed line and the dotted lines 
are, respectively, the best-fit observed relation (equation 1331 
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with logc C 2 = 0.54 and 8 = 1.03), and the associated scatter 
in 8 and logc e 2- In Fig. [9] we report plots for the SIS, LTM 
and acNFW+S models. We note that the acNFW+S models 
considered in this case have fixed value of /„ (/„ = 0.1 or 
/* = 0.02), but the value of £ is a function of m: for given 
m we obtain R D from the observed relation (|3ip and then £ 
from equation p4)l . fixing C200 = 8. 

From the diagrams in Fig. [9] it is apparent that the 
SIS models behave very differently from the LTM and ac- 
NFW+S models (NFW+S models, which are not plotted, 
behave similarly to acNFW models). Let us focus first on 
the case with r a /7? e independent of m (left-hand and cen- 
tral panels in Fig.©. For the SIS model the ratio M^ 2 /M^ m 
gradually increases with M^ 2 , while for all other models the 
ratio M^ 2 /M^ 2 m increases with M^ 2 at low masses and de- 
creases with M£ 2 at high masses (in sharp contrast with the 
tilt predicted for SIS models), and the variations with M£ 2 
are stronger than in the case of the SIS models. Quantita- 
tively, at high masses, the predicted slope 8 of equation (|33[) 
is 8 > 1 for SIS models and 8 < 1 for the other models. 
When r a /R c depends on m (right-hand panels in Fig. [5]), 
the dynamical non-homology introduces additional tilt, in 
the sense that the predicted value of 8 at high masses be- 
comes smaller, giving 8 < 1 also for SIS models. However, 
also in this case significantly less tilt is predicted for the SIS 
models with respect to the other models. 

In conclusion, the effect of structural non-homology on 
the tilt of the MP is of order of a tenth of a dex in the mass 
ratio and thus measurable if one had a sample comparable 
in size and quality to SLACS, covering a further decade 
down in galaxy masses. Interestingly, the tilt of the MP is 
measurably different depending on whether or not the total 
mass profile is well represented by a SIS. The tilt of the 
MP appears thus to be a powerful diagnostic of the internal 
structure of early-type galaxies. 



5 SUMMARY AND CONCLUSIONS 

In the range a C 2 = 175 — 400 km s" 1 the MP of early-type 
galaxies has no significant tilt and small associated scat- 
ter. This means that the dimensionless structure parameter 
c e 2 defined in equation (0 is a nearly universal constant. 
In other words, the range of values of c C 2 "allowed" by the 
observational data is remarkably small. Even for spherical 
galaxy models, c e 2 is expected to depend on the stellar den- 
sity profile, orbital anisotropy of stars, and total (dark plus 
luminous) mass distribution. In this paper we explored the 
constraints posed by the existence of the MP on several rel- 
evant families of galaxy models(f| 

Limiting to spherically symmetric models, we found 
that c e 2 is not very strongly dependent on galaxy struc- 
ture and kinematics, so a relatively wide class of models 
have values of c G 2 within the observed range. Therefore, 
strictly speaking, the massive early-type galaxies lying on 
the MP are not necessarily structurally and dynamical ho- 
mologous. However, not all the studied models behave in 

4 Throughout the present paper we considered the MP in the 
standard contex t of N ewtonian gravity with dark matter. See 
ISanders fc L and (2008j) for an interpretation of the MP in the 
context of Modified Newtonian Dynamics. 



the same way when compared to the observational data. 
Models in which the total density profile is a SIS are con- 
sistent with the observed range of c C 2 for a wide class of 
stellar density profiles, and only models with extremely ra- 
dial or tangential anisotropies are excluded. The light-traces- 
mass hypothesis is not excluded by the observational con- 
straints here considered, apart for the case of high-m Sersic 
models, which cannot be reconciled with the MP within 
the observed scatter. (However, LTM models are known 
to fail other observational constraints: see Section f2.3p . We 
also considered cosmologically-motivated models with NFW 
dark- matter halos (with or without adiabatic compression), 
finding that they are consistent with the MP only for a rel- 
atively limited range of values of their parameters, so a de- 
gree of fine-tuning between light profile and anisotropy is 
required. Among these NFW models, those with adiabati- 
cally contracted halos and those that are baryon dominated 
seem to require slightly less fine tuning than those with non- 
contracted halos and those that are dark-matter dominated. 

This work has focused on the average value of c C 2 in 
the SLACS sample, along with its intrinsic scatter. With 
the exception of Section [4] we have not explored the impli- 
cations of the fact that this intrinsic scatter is not correlated 
with either mass or with the ratio of Einstein radius to R e 
ijBolton et al.ll2008bh . These observational results indicate a 
degree of structural homogeneity across a range in mass. In 
future works, we will explore these mass-dependent results 
in the context of mass-dynamical models such as those con- 
sidered here. We also plan to refine these analyses based on 
the results of forthcoming velocity-dispersion measurements 
of higher signal-to-noise ratio and in smaller and more uni- 
form spatial apertures. 

We also explored the possibility that the observed de- 
pendence of the Sersic index m on the galaxy luminosity 
could tilt the MP when a sufficiently large mass range is 
considered. In this respect, SIS models behave differently 
from all other models: a slightly tilted MP is predicted in 
the early type galaxies have SIS total density distribution, 
while a "bent" MP is predicted in all the other explored 
cases. The effect is large enough to be measurable with sam- 
ple of lenses comparable to SLACS in size and quality and 
extending a further decade in galaxy mass. 

In conclusion, our results are consistent with the hy- 
pothesis that massive early-type galaxies have isothermal 
(oc r~ 2 ) total mass density distribution, though alternative 
hypotheses cannot be excluded on the basis of the existence 
of the MP alone, although in some cases they require a de- 
gree of fine tuning. In any case, the process of formation 
of early-type galaxies lead to systems with a combination 
of total mass distribution, luminosity profile, and orbital 
anisotropy such that they lie close to the MP. It will be 
interesting to quantify whether the observed fine tuning is 
quantitatively consistent with the range of simulated proper- 
ties of early-type galaxies in the standard hierarchical model 
of galaxy formation. 
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Figure 3. Contours of constant logc c 2 in the planes /3~7, f3-m (left, for /3 = const models), and log(r a /iJ )-7, log(r a /-Re)-"i (right, 
for OM models), for systems with singular isothermal sphere total mass distribution. Dark-shaded region s correspond to non-consistent 
models. Light-shaded regions are outside the observed range (logc e 2> = 0.53 ± 0.057 teolton et al.ll2008bh . 
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Figure 4. Same as Fig. [3] but for light-traces-mass models. 
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Figure 5. Same as Fig. [3j but for NFW plus stars models with /* = 0.02. 
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Figure 6. Same as Fig. [3] but for NFW plus stars models with /* = 0.1. 
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Figure 7. Same as Fig. |3] but for adiabatically contracted NFW plus stars models with /„ = 0.02. 
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Figure 8. Same as Fig. [3] but for adiabatically contracted NFW plus stars models with /* = 0.1. 
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Figure 9. Ratio between the total projected mass within iJ c /2 (Mf 2 ) and the dimensional mass M c d 2 im = R e al 2 /(2G) as a function of 
M p 2 for R}/ m Sersic spherical OM galaxy models with different total mass distributio ns and values of r a /i? e . In each panel the dashed 
line is the best-fit observed relation with the associated scatter (dotted lines) found bv lBolton et alj l l2008bl , see text for details). 



